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Abstract 

Electromagnetic properties of quark-like particles are examined in a classical field model involving 
extended dual electromagnetic fields. These can have fractional charges and a confining potential 
that derives essentially completely from a short-range weaker potential. The combined potentials 
exhibit an asymptotically free spherical surface and contribute to the masses of the particles. The 
quarks are shown to have an intrinsic symmetry that describes their structures in hadrons. Multi- 
quark solutions are easily obtained for both stable and unstable particles. Each quark can undergo 
simple harmonic motion in a range of frequencies. 

PACS numbers: 12.39.-x 



1 Introduction 

In this note a classical electromagnetic model of particles is introduced that have a number of quark-like 
properties. Starting with the Lienard-Wiechert exact solutions of Maxwell's equationj^ of a point charge 
q in arbitrary motion, it can be shown that the following relation between the fields and potentials holds 
everywhere: [e = c = 1] 

{B^ - E'') = -q-\A^ - <f'f (1.1) 

With F^^ = dfj_A,j — 5yj4^, 3-potential vector A, ^4 = i4> and the relation F^,^F^i, = 2{B'^ — E'^), (1.1) 
can be written 

F^,,F^, = -2q-''{A^A^f (1.2) 
The Lienard-Wiechert solutions also satisfy the orthogonal relation E • B = everywhere, i.e., 

*F^uF^, = (1.3) 



*Mailing address: 304-2323 Hamiota St. Victoria, BC. CA VSR 2N1; email: hschiffOshaw.ca. 

^See, for example, J. D. Jackson, Classical Electrodynamics, (John Wiley & Sons, Inc. 1975, 2nd ed) p. 654 ff. 
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In (1.3) the familiar notation = 1/2 ie^napFa^ is used for the dual of F^i, using the Levi-Civita 
tensor e^j,i,ap- Due to the opposite symmetries of Ff^i, and A^A^,, (1.2) can be written, 



Defining the mixed tensor 
(1.4) can also be written 



{F^, + ^Iq-^A^A.f = (1.4) 
W^, = F^, + ^2q-^Ai,A,, (1.5) 



VFm-^m- = (1.6) 
Further, noting that W^,/ = since Sni/a^AaAjj = 0, (1.3) can be written, 

*W^,,W^, = 2E ■ B = (1.7) 

The expressions (1.6) and (1.7) show that W^i, (1.5) satisfies formally the conditions for a null field 
(although the term 'null field' applies strictly to an anti-symmetric tensor we have taken the liberty to 
use it here for the mixed tensor) that holds everywhere in space, as compared for example to solutions 
Fjj^i, to Maxwell's equations where, in general, only the radiation field is a null field. This interesting 
property of W^j, a combination of the exact solutions for the electric and magnetic fields and potentials 
of Maxwell's equation for a point charge, arouses curiosity as to whether the expression M^i_ii, may exhibit 
other features of interest. To this end we consider the non-gauge invariant equation d^W^i, = — 47r5i/, 

^^,iF^u + ^2q-^A^A,) = -4n% (1.8) 

where % is an external localized 4-current density {^3, if)) that can be chosen freely. The total cur- 
rent density thus consists of the localized external current % density plus the self-field contribution 
^2{A7rq)-^d^{A^A,). 

2 Solutions of (1.8) 

We will examine static, radially symmetric solutions of (1.8), 

y/2q-^V ■ {AA) = -ATTV{r) (2.1) 

V • E - ^J2q-^V ■ ((^A) = AirXr) (2.2) 

The asymptotic solution for (2.1), 

V • (AA) = 0, (2.3) 

A = a/r (a arbitrary) (2-4) 
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Defining 

7 = ^1q-^a, (2.5) 
A = 7g/V2r (2.6) 
Asymptotically, using (2.5), (2.2) can be written, 

^ + 2:^ = (2.7) 
For any value of 7, except 7 = 1, treated separately below, there are two indicial solutions of (2.7), 

</>! = h/r (2.8) 

(t>2 = cr-^ (2.9) 

A Coulomb potential and a possible absolute confining potential for 7 < 0, where c may be positive 
or negative. For any value of 7 the total charge of the external source in (2.2), using the divergence 
theorem, 

J'pd^x = {l-j)b (2.10) 

For 7 = 1 the two indicial solutions of (2.7) merge to a single Coulomb potential so a second solution 
is needed. This is easily seen to be given by rcf) ~ logr. We note that for 7 = 1 the total external charge 
is zero and one can choose its charge density to be zero. These solutions however will not be used here. 
In the following the particular solution to (2.2) for which (2.8) is the asymptotic solution will be called 
(j)p. 

Using the divergence theorem in (2.1) the total 3-current of the source is 

j V{r)(fx = -l/47r j ^2q-^A\'^dn = -7^g/V2 (2.11) 

3 Null Condition and Fractional Charges 

As A and go as 1/r asymptotically, the application of (1.2) is suggested to these solutions, leading to 
an equation for the charge b. Thus, with B = 0, (1.2) becomes 

-E2 + g-2(A2- 02)2 = (3.1) 

For each sign of q one obtains a quadratic equation for the Coulomb charge b in (2.8), 

b^±qb- -i^ijl = (3.2) 
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Solutions 6_, with the minus sign in (3.2), have opposite signs to those with the plus sign 6+, equivalent 
to g — )• —q. Thus for the quadratic with the minus sign, 

b'^ -qb- = 0, (3.3) 

b- has the two solutions, 

6_(±) = g[l±V(l + 27')]/2 (3.4) 
In (3.4) choosing 7^ = 4 and g = ±1/3, one gets 

6_(±) = ±(2/3, -1/3), b+ = -b_ (3.5) 

(With b = ±1/3, ±2/3 and 7 = —2 (negative for confinement), the external charge (2.10) is ±1,±2 
respectively. Other possible fractional values for b can be obtained, for example with 7^ = 4 and 
q = ±2/3, 6_ = ±(4/3, -2/3), 6+ = -6_). 

One may also note that with 7^ = 4, choosing only the minus sign on the right hand side of (3.4) 
with g = ±1 one gets the two values b = —1 and ±1. 

4 Lagrangian 

A Lagrangian for (1.8) cannot be obtained using only the field components in that equation. The 
nonlinear addition to the second order Maxwell terms of the potentials is in first order, a problem 
similar to that of Newton's equations of motion with a first order dissipative term. This is resolved for 
(1.8) with the addition of an electromagnetic fielc]^ G^jy and its associated potential VJ^, 

G^. = d^V, - d,V^ (4.1) 

The following Lagrangian includes F^y,Ay and G^y,Vy with their respective currents '^ju and ^u- 

L = -l/87r[l/2G^,F^, ± ^2q-^V,d^{A^A,) + AW,% + AttA,%] (4.2) 

Variation of L with respect to gives (1.8), while variation with respect to A^ yields a linear equation 
for V^. 

d^G^, - ^2q~^A^{d^V, + d,V^) = -A7T% (4.3) 
The external source may be chosen freely. The total current density 

= -^2{A7:q)-^A^{d^V, + d,V^) + % (4.4) 

^See, for example, Morse & Feshbach, Methods of Theoretical Physics (McGraw-Hill Book Company, Inc., 1953), Vol. 
1, p. 298, 313. Where they discuss introducing mirror image systems for dissipative systems. 
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is conserved, and for radial solutions considered here the total 3-current density J = 0. 
Prom (4.3) the equation for V4 = iil^, becomes 

Both particular and homogeneous solutions of (4.5) will be considered. Before dealing with general 
solutions to (2.2) and (4.5) in Section 6, it is useful to look first at an explicit solution for a single- 
quark. First, for the homogeneous function V'h of (4.5), replacing ■\/2q~^A by A*, 

^ + ^-A*#^ = (4.6) 
ar^ rdr ar 

Using the asymptotic value of A in (2.6) with 7 = —2, (4.6) becomes, 

-7^ + — = (4.7 

ctr^ rar 

It follows that for all asymptotic tp^ solutions, 

V'/j ~ + const. (4.8) 
To find A* in (4.6) we choose for i/jh, to be defined as ipx, the simple monotonic function 

iPh = i^x = AsV(s + arf + XsH (4.9) 

A is a charge plus an amplitude, s is a scale parameter and a a multiplier for separating distance between 
s and r. An arbitrary constant t has been added to tpx with dimension (One expects s to be 

small compared to the position of the quark, i.e., a <^1). Using (4.9) one finds for the 3-potential A* 
in (4.6), 

A* = 2/r-Aa/{s + ar) (4.10) 

With (4.10) in (2.2), which is a complete divergence, putting the divergence constant K in K/r"^ equal 
to zero, one obtains the following homogeneous solution for the confining potential, 

(j)k = k{s + arf/s^r'^ (4.11) 

A; is a charge plus an amplitude. (j)k corresponds asymptotically to (2.9) {k/s^ replaces c) and diverges 
as 1/r^ at the origin with a minimum value of 16ka^/s at ar = s. 

Substituting (4.10) in (2.1) one finds for the external 3-current density 

J{r) = -2,/2q[S{r) - a^s(s - ar)/7rr^(s + arf] (4.12) 

Spatial integration of (4.12) is in accordance with (2.11). 
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For the 3-potential V in (4.3) one can choose the external 3-current density *J = 0, the equation for 
the V would then become, 

2A*^ = ^ (4.13) 
dr dr 

However, since V will have no bearing on the physical quantities we consider, its solution will not be of 
interest. For the particular solution (j)p of (2.2) we are free to choose, 

= 6aV/(s + ar)^ (4.14) 

For energy considerations later the particular solution ^/^p of (4.5) is taken to be the same function as 
(4.14), 

iljp = (j)p = ba^r'^/{s + arf (4.15) 

In (4.5), as ^pp, (pp and ^* go as 1/r asymptotically, both the field and external charge densities are seen 
to have a 1/r^ singularity with opposite signs that cancel for a non-singular total charge density. 

5 Energy and Mass 

From the Lagrangian (4.2) the canonical stress-energy tensor is 

dL dL 

To represent the energy of a particle, we choose for simplicity, the stress-energy tensor of the elec- 
tromagnetic fields of (5.1), Fafj, and G^^, 

Qa^ = l/8n[Ga^,Fp^ + Fa^G0^ - l/2{G^,G^,)6ap] (5.2) 

Further, in (5.2) only the particular solutions ijjp and <pp are chosen for Ga^ and Fq,^ and as these 
solutions are now equal, = F^jj,. This leads to a familiar expression for this energy density, 

- G44 = l/87r[-2F4^F4^ + 1/2{F^,F^,] = l/SnidcPp/drf (5.3) 

with integrated Coulomb energy 



Co 



1/2 j {d(l)p/drfr'^dr (5.4) 



For the contributions to the energy by the homogeneous potentials (pk and ipx we choose these 
energies to be the interactions of a quark charge h with the potentials, hcpk and bip\. Thus, when added 
to (5.4), the total energy for the single-quark, (4.9) and (4.11), becomes, 

^ = ^^ + ^j^ + ^y^ = 0.043a6Vs + bk{s + ar)^/sV + h\s^/{s + arf + h\sH (5.5) 
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For the numerical value in (5.5) the function (f)p in (4.14) was used. For the mass one needs the minimum 
of + plus the Coulomb energy. In general, the result for the minimum will depend on k, A, and 
a and will be different from the minimum — such as IGbka^ /s at ar = s for the single quark. The 
spherical surface corresponding to the minimum energy is evidently one of 'asymptotic freedom'. The 
energies and masses for the single-quark and multi-quarks are dealt with in the next section. 



6 Multi-Quark Solutions 

General solutions for (f)^ in (2.2) requires A* from (4.6), 

^* = V'AM + 2/r, (6.1 
The dimensional constants Xs^ cancel in i/j'^/i/j'^, replacing ipxjXs^ by the capital one can write 

V'AM = *AM = ^^log(|^l|)/dr 

and 



A* = dlog{r'^\^'x\)/dr 
In (2.2), with the equation for (f)]^ a complete divergence, replacing (t>k/{k/s'^) by ^k-, 

d^k/dr + = Ks^/r'^ 

In (6.4), has been added for a dimensionless K. With A* from (6.3), (6.4) becomes 

d^k/dr + ^k[d{r^\'^'x\)/dr]/r^\'i>'x\ = Ks^/r^ 



Multiplying by r"'!*';^!. 



r^\^'x\d^k/dr + $fcd(r^|*';,|)/dr = Ks^'l^'x 



2hT,/ 



„3hT,/ 



d{<l^krVx\)/dr = Ks'^\^'x 



^k = Ks-" 



\^'x\dr 



/r'^\^'x\ + C/r'\^!'x 



2|,T,/ 



With 



j \^'x\dr /r^l^Al = j ^'xdr 



(6.2) 
(6.3) 



(6.4) 

(6.5) 

(6.6) 
(6.7) 
(6.8) 

(6.9) 
(6.10) 



= Ks^i^x + li)/r''^'x + CA'I^aI 

The constant /3 (dimension 1/s^), consists of the integration constant plus the arbitrary constant 
added to ^x- Putting Ks^P = P, one obtains the three distinct terms. 



= Ks^^x/r^'^'x + P/r'^'^'x + C/r'^\% 



3.11) 



We note that the K and P terms change sign in crossing an extremum of ^x. With different various 
multi-quark solutions follow from this general solution to the confining potential, defined by the relative 
signs and magnitudes of the three dimensionless parameters K, P and C. Some examples are shown 
below. 

From (6.3), asymptotically A* —2/r and near r = 0,A* ^ {n + l)/r for ^ r'^,n> 1, between 
these points A* diverges at the extrema of ^x- Near the origin in (6.4), with The K term from (6.11), 
one finds readily the necessary equality with K/r"^; including all term in (6.11), goes as ztl/r""^^ 
and diverges at the extrema of ^x- Asymptotically the K term goes as ±l/r and the P and C terms 
go as ibr^. We will consider a few choices for the magnitudes and signs of the three parameters K,P 
and C, for a variety of hadrons. To summarize, for all solutions of (6.11) and (6.3): 

(a) $fc ^ ibr^ asymptotically 

(b) ~ ibl/r"'''"^ near the origin, n > 1 

(c) <I>fc and A* diverge at the extrema of \ (6-12) 

(d) A* —2/r asymptotically and as (n + l)/r near the origin 

The behavior in (6.12) is readily seen to apply to the single-quark solutions (4.9), (4.10) and (4.11), 
except for (c), as (4.9) has no maximum or minimum. 
For the simplest case in (6.11), where K = 0, 

= C/r^\^'x\ (6.13) 

From (6.3) and (6.13), 

A* = -(l/$fe)d$fe/dr (6.14) 

In this case A* = aX the extrema of 

For a 2-quark system a function ipx is needed with one maximum. The 0^ divergence at this 
maximum (6.12c) divides space into two regions in each of which (p^ is contained by two divergences, 
between the origin and the central divergence in the first region and between the central divergence and 
asymptotically in the second region. As an example of a simple ipx with one maximum, applied to a 
charged 2-quark system, we consider 

'iIjx = \s^r/{s + arf (6.15) 

^'x = Xs'^{s-3ar)/{s + arf (6.16) 

(j)k = k/s^{s + arf/r'^\s-3ar\ (6.17) 

In (6.17), (pk was obtained using (6.13), with C included with the amplitude in k. At the maximum 
of dtpx/df, where r = s/3a, the (pk divergence there divides the space into the two regions, r < s/3a 
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and r > s/Sa. In region(l), (pk has a divergence at r = and at r = s/Sa, while in region(2), (f)k has a 
divergence at r = s/3a, and asymptotically. The zero derivative of (6.17) results in a quadratic for the 
two (f)k minima, r^m = 0.18s/a and 1.82s /a, straddling the divergence at r = s/Sa. In each region the 
quark mass, with appropriate b, obtains from the minimum of the energy 6(0^ + 

In the region of a quark in a multi-quark system, dtp^/dr is either positive or negative, so the energy 
minimum there, from d[b{(f)k + il)\)\/dr = 0, can be written, 

d(j)k/dr = -di^x/dr = ±\^'x\, (6.18) 

where the + or — signs refer to ■^^ < or > respectively. With K small enough to be ignored, 
introducing the dimensional factors and using (6.13), (6.18) can be written, 

d^k/dr = ±s^X\'^'^\/k = ±s^X/kr^^k (6.19) 

^kd^k/dr = ±s^X/kr^ (6.20) 

(To illustrate, we take A; > 0, A > 0). An r-valuc, r„j, that satisfies (6.20) minimizes the energy 
Cfc + Ca between two divergences. When ^'^ < and d^k/dr > 0, the 1/r^ curve in (6.20) intersects the 
right half of $fc and the actual minimum, r„i, for that quark moves to the right of the minimum, 
''kmin- (For tlic singlc-quark solution above whose minimum is at r = s/a and < 0, the actual 
minimum > s/a). Succinctly: 

< ^ d^k/dr >0^rm> rkmin (6.21) 

> ^ d^k/dr <0^rm< rkmin (6.22) 

To apply the use of (6.20) for determining the minimum energy (for simplicity) we use the single- 
quark solution (4.11). With = ns/a one obtains the following equation for n, ((6.21) applies), 

2(n + l)^(n - l)/n^ = +X/ka^ (6.23) 

Choosing, for example, the simple values X/k = 10"^ say and a = 10~^, X/ka^ = 10^, the solution is 
n = 12.6. We would want to be less than one fermi, 10~^^m, so with = ns/a < 10~^^,s < 
W^^a/n, thus s < 10~-'^^m. (For X/k = 1 and a = 10""^, n is again 12.6 giving s < 10~^^m). 

A neutral system would require both P and C terms in (6.11). With small enough K there are two 
possible choices for P and C where the first > and the second < 0, or vice versa: 

[P>0,C <0:P> -C], or [P < 0, C > : -P > C] (6.24) 

(For positive or negative charged systems appropriate choices for P and C are easily made). 
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For a 3-quark solution, with the 3rd ^'^ < 0, '^x must have one minimum followed by one maximum. 
These two extrema produce two divergences in the central region and thus three extrema in for three 
quarks. In region(l) < 0, in region(2) > and in region(3) < 0. For a proton, say, two positive 
and one negative 0^ are needed. For a small enough K (6.11) and A; > 0, with P < 0, C > 0, —P > C, 
the first 0^ is positive, the second is negative and the third is positive, so the energy, 6^^, of each quark 
in its region is positive. Another choice is P < 0, C < and —P > —C. The two negative slopes 
represent an up quarks and the one positive slope a down quark, (the d quark always lies be between 
the u quarks). Formally, n extrema result in n + 1 quarks with an infinite barrier between adjacent 
quarks, arranged in 'shells' of asymptotically free concentric spherical surfaces, one quark per 'shell'. A 
possible symmetry that groups the quarks into hadrons is discussed in the Symmetry section. 

In the above examples, in each case, the system is completely confined. With the inclusion of a large 
enough K term however an unstable quark can result. Choosing with a negative slope for the second 
quark of a 2-quark particle and choosing C > with P > C, (the P term changes sign compared to its 
sign in the first quark) the sum of both P and C terms in (6.11) produce a negative divergence at the 
^\ maximum and another negative divergence asymptotically. Including now the K term contribution, 
the second quark's in (6.11) can be written (the K term also changes sign), 

$fe = -Ks=^*aA^|^aI -{P- C)A^I^aI (6-25) 

With K < and — -f^s^^'Amax > P — C, the total divergence at ^'Amax will be positive for a 
positively charged particle, and asymptotically where the K term ^ as +l/r, the divergence of 
will be negative, eliminating absolute confinement. Shown below, solutions exist for the second 
quark that have a minimum followed by a maximum and then the — divergence, giving only partial 
confinement. Using (6.25), with a charge b for the second quark, the energy minimum (and maximum) 
is obtained from the zero derivative of the interaction energy, 

^k + ^\ = K^'k + V'A + const) 

= bk/s^[-KsHx/r'^\^'x\ - {P - C)A^I^aI] + bXs^'^x + bXs^t (6-26) 

where (P — C) > 0. (For the first quark the energy will have +(P + C) in (6.26)). To show solutions 
for an unstable quark we choose *a from the 2-quark function (6.15). Dividing (6.26) by hk/s, putting 
X = ar/s, the energy can be expressed in the dimensionless form, 

E = D{1 + x)/x{?,x - 1) - U{1 + xf/x^{?,x - 1) + Nx/{1 + xf + T (6.27) 

where D = —aK,U = a'^{P — C),N = X/ak and T = s^Xt/k. (At this second quark, the U term 
confirms that asymptotically, E diverges as — x^). For a positive divergence at a; = 1/3 (the maximum 
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of from (6.16)) D/U must be greater than 9.48. Without loss of generahty we put U = 1. The zero 
derivative of (6.27) gives the following relation between D and A?^ for x > 1/3, 

^ 2(l + x)4(3x^-6x + l) (3x-l)W 

x(-3x2-6x + l) (l + x)5(-3x2-6x + l) ^ ' 

For a few solutions, Table 1, we chose, arbitrarily, a specific value for D and a specific value for x to 
be an extremum; the following values were obtained for A?", Xmin, ^Cmax, -Emin and -E^max- The variation of 

Table 1: Unstable Solutions 

D N Xjnin ^max -E^min -E^max (-E^max -E^min) 



10 


-83.4 


0.8 


3/2 


-12.3 


-10.8 


1.5 


10 


-152 


0.7 


2.0 


-17 


-12.9 


4.1 


10 


-488 


0.5 


3.0 


-49 


-18.3 


30.7 


102 


-1163 


2.0 


2.5 


-10.93 


-10.77 


0.16 


10^ 


-10390 


2.3 


3.0 


29.1 


30.7 


1.6 



depths of the various minima, (.Emax — ^min), range from about 0.16 to 31, reflecting various lifetimes 
for the different particles. (The initial chosen extremum turned out to be a maximum except the one 
for D = 102). To each set with negative energy at the minimum a value of T must be added to shift 
to a positive energy. There is no solution with 7 = — 2 when C = P{U = 0), for then (f)^ would go as 
1/r asymptotically. With the decay of an unstable charged 2-quark particle, 0^ — >■ 0, leaving (f) = (f)p 
and = (f)p + Without confinement the interaction energies b(f)k and bipx are superseded and the 
emitted free particle, with unit charge and corresponding equal (f)p and ipp (a possible leptonic decay), 
would have a mass obtained from (5.2), 

m = 1/2 J [{d^p/drf + {d(f)p/dr){d'ipx/dr)]r^dr, (6.29) 
showing a contribution to the mass from the weak potential ipx. 



7 Simple Harmonic Motion 

It is instructive to examine a quark's simple harmonic motion near the bottom of its energy-well, its 
asymptotically free spherical surface. The force on a quark, —d^/dr, for a small displacement e from 
Vm at the energy minimum ^(r^) is 

F = ~e{rm + e) = -eirm)e + 0{e^) + ... (7.1) 
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With mass m = ^{rm)/(?, (Coulomb energy is assumed to be small), the equation 

m(fe/dt^ = F 



(7.2) 



becomes, 

Se/dt" = -[^"{rm)/i{rm)]c^e (7.3) 

Hence, the angular frequency 

= r{r^)/ar^)\^'^c (7.4) 
Again, using the single-quark solution (4.11), one finds for the new minimum at r = ns/a, 

= 2hka^/s\n + 1)^(^2 - 2n + 2,)a^ + &n^\/k]/n^{n + 1)^ (7.5) 

^ = hk/s[{n + l^a^ + n^\/k]/n^{n + 1)^ (7.6) 

^ 2a^[(n + l)^(n^ - 2n + 3)a^ + Gn^A/fc] 
^ s2[(n + l)7a2 + n2A/fe]n2(n + l)2 ^ 

Using the values from the example for (6.23) where \/k = a = 10~^ and n = 12.6, one finds the 
frequency from (7.7) and (7.4), 

uj ~ O.lca/s = 3 X 10^/s > 3 x lO^^/sec, (7.8) 

corresponding to a photonic energy huj/2Tr w 0.2GeV. (For A/A: = 1 and a = 10~2 where n is again 
12.6,/iw/27r fsi 0.02GeF). The and Va energies from (7.6) with n = 12.6 and A/A; = lO'^ plus the 
Coulomb energy are: 

2 X 10-^6A/s, .^A ~ 4 X 10-^6A/s = 2 x lO'^^fc, ^coui ~ 4 x (7.9) 

Assuming k > b, the energy is dominated by ^k- Writing b = b'e,k = k'e (with k' the amplitude, h' 
either 1/3 or 2/3 and e the electronic charge) and s = 10~^^m,^fe ^ b'k'3MeV. 



8 Symmetry 

These quarks consist of products of two equal particular solutions, (^p and tpp, (shown symbolically 
distinct for clarity — subscript p will be suppressed). With the two different q values ±1/3 for which the 
same fractional charges obtain (3.4), (3.5) four different (pip products can be made, such as (p{l/3)'ip{—l/3) 
etc., representing 4 possible quarks (or anti-quarks) with the same charge. For brevity we show an up 
quark, say, {(A2/3(l/3)V'2/3(~l/3)} with the g-pair (1/3,-1/3) as u(+,— ) similarly for other quarks. 
The pairing of the q values with the (ptp combination implies that there are effectively 4 distinct g-pairs 
(+, +), (— , — ), (+, — ), (— , +) with the obvious property that the sum of their 8 individual q^s is zero; 
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but as we see next, zero sum q''s (g-neutrality) is a necessary and sufficient condition for tfie formation 
of baryons and mesons (evoking QCD's color neutrality). 

From the 4 g-pairs a baryon's three quarks can always have three g-pairs with zero sum g's consisting 
of the two g-pairs (+,+),(—,—) plus, randomly, either the g-pair (+,— ) or (—,+). A proton, for 
example, could change continually between [u{—,—),u{+,—),d{+,+)] and [u{—,+),u{+,+),d{—,—)] 
(order of g-pairs arbitrary). Mesons can also use all 4 g-pairs in their sets of [(quark), (anti-quark)] 
combinations, [(+, +), (— , — )] and [(+, — ), (— , +)] plus the two sets with the signs reversed, all sets with 
zero sum g's. (How g-pairs may change is discussed below in Conclusions). 

With the bi-functional form of quarks as the basis for the effective non-commutativity of g-pairs 
(+, — ) and (— , +), the 4 distinct g-pairs (+, +), (— , — ), (+, — ), (— , +) show that g-neutrality can occur 
in 4, 3, 2 and 1 g-pairs, with possible relevance to multi-quark particles beyond hadrons as well as 
single-quark particles; though the latter are consistent with the single-quark solutions obtained above, 
these have not been observed [1]. 

9 Conclusions 

The Lagrangian (4.2), although missing initially the usual kinetic terms F^^Ff^u and G^yG^y, was 
restored to the degree that the choice of equal particular solutions of the two field equations merged the 
corresponding F^^ and G^i^. This led to an electromagnetic stress-energy tensor with an energy density 
corresponding to the Coulomb energy density of a single charge (5.3) — raising inevitably, the intriguing 
question of a possible quark structure. 

The results, such as fractional charges and absolute confinement, including asymptotic freedom, are 
in accord with those of QCD; as well, the manifest intrinsic symmetry adds support to the model. The 
g-pairs, combined with the bifunctional form of the quarks establish the quarks' symmetry, but are not 
involved in any direct interaction. 

Solutions with multi-quarks are easily obtained describing both stable and unstable particles. The 
significant result where the strong confining potential 0^ is completely determined by the short-range 
weak potential highlights the unity between them and electromagnetic potentials, emphasizing also 
the primacy of the weak potential with respect to the derivative strong confining potential. Prom the 
bottom of the energy-well on its spherical surface a quark can move radially in simple harmonic motion. 
This of course does not preclude a possible additional force-free motion on a spherical surface. 

How the g-pairs in hadrons can continue to change is not, in general, within the purview of these static 
solutions. However, with expected time dependence, the q values in (1.8) could be viewed as changing 
rapidly and continually between 1/3 and —1/3, leading to random changes in the g-pairs. Classically, 
this 'virtual' process should not be observable; the average q value would be zero and (1.8) would revert 
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to the familiar Maxwell equations. However, because of the success of the static solutions in describing 
important quark properties, we suggest that they could be interpreted as displaying 'snapshots' of this 
process. 

Overall, this presents a different and possibly fruitful perspective on quarks within a framework of 
classical fields where all forces and interactions are electromagnetic. 
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